Critical Currents and Vortex States at Fractional Matching Fields in Superconductors 

with Periodic Pinning 
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We study vortex states and dynamics in 2D superconductors witli periodic pinning at fractional sub- 
matching fields using numerical simulations. For square pinning arrays we show that ordered vortex 
states form at 1/1,1/2, and 1/4 filling fractions while only partially ordered states form at other 
filling fractions, such as 1/3 and 1/5, in agreement with recent imaging experiments. For triangular 
pinning arrays we observe matching effects at filling fractions of 1/1, 6/7, 2/3, 1/3, 1/4, 1/6, and 1/7. 
For both square and triangular pinning arrays we also find that, for certain sub-matching fillings, 
vortex configurations depend on pinning strength. For weak pinning, ordering in which a portion 
of the vortices are positioned between pinning sites can occur. Depinning of the vortices at the 
matching fields, where the vortices are ordered, is elastic while at the incommensurate fields the 
motion is plastic. At the incommensurate fields, as the applied driving force is increased, there can 
be a transition to elastic flow where the vortices move along the pinning sites in ID channels and 
a reordering transition to a triangular or distorted triangular lattice. We also discuss the current- 
voltage curves and how they relate to the vortex ordering at commensurate and incommensurate 
fields. 



I. INTRODUCTION 

Vortex lattices interacting with quenched disorder rep- 
resent an ideal system in which to study various static 
and dynamic phases of elastic lattices in disordered me- 
dia, which can be found in a large class of condensed 
matter systems. Pinning of vortices is also of practical 
importance since superconductors are required to main- 
tain high critical currents for various potential technolog- 
ical applications. There have been considerable efforts to 
enhance the pinning properties of vortices in supercon- 
ductors by introducing point, columnar, and splay defects 
using electron, proton and heavy ion irradiation j^]. A 
particularly promising technique has been to lithograph- 
ically create nano-structured periodic defect arrays. Ex- 
periments with arrays of holes ^-p^] and magnetic dots 
[14 |l9|, as well as theory and numerical simulations 
[21 1, have shown that a variety of commensurability ef- 
fects occur when the density of vortices equals an integer 
multiple of the density of pinning sites such that a well 
ordered vortex crystal can be stabilized and pinned. Sev- 
eral experiments have provided evidence that for small 
pinning sites, beyond the first matching field, a portion 
of the vortices will sit at the interstitial regions between 
the pinning sites ||6|,[l0|. 

Besides matching effects at integer filling frac- 
tions, commensurability effects at non-integer and sub- 
matching fields have also been observed directly with 
Lorentz microscopy and indirectly with magnetiza- 
tion and transport measurements . Experiments by 



Metlushko et al. observed sub-matching effects down to 
1/16 filling fractions ||]. These same experiments did not 
show any particular fractional matching for B / Bfj, > 1.0. 
In other experiments, fractional matching has been seen 
at filling fractions of 3/2 and 5/2 |l|]. These fractional 
matching effects are generally less pronounced than those 
of the integer matching fields. 

Recent scanning-Hall probe measurements have also 
directly imaged vortex configurations at various sub- 
matching fields ||l^,|l3|. These measurements find that 
some configurations are only partially ordered while oth- 
ers, including the integer matching fields, are almost com- 
pletely ordered. Experiments by Field et al., who image 
up to 5,000 vortices, have found that at the 1/2 matching 
field very large regions of ordered vortex checkerboard 
states can be observed with well defined domain walls 
separating different orientations At other fields, 

such as 1/3, only domains of diagonally ordered regions 
are observed, while at 1/4, 2/5, and 1/5 even smaller 
regions of order are observed. 

Previous simulations of vortices interacting with pe- 
riodic pinning arrays have focused on vortex matching 
at integer matching fields and have found that differ- 
ent types of vortex crystals can be stabilized at dif- 
ferent matching fields in agreement with experiments 
[^ . The case of vortex matching at fractional fields for 
square and triangular pinning arrays has not been ad- 
dressed. Although commensurate vortex configurations 
in Josephson-junction arrays and repulsive particles 
on lattices at fractional fillings have been studied, 
these models do not allow the vortices to sit at intersti- 
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tial regions between the pinning sites as is possible in 
experiments with periodic pinning, magnetic dot, or hole 
arrays. For B/B^ > 1 it is clear that interstitial pin- 
ning of vortices is possible. In this work we show that 
for certain pinning strengths and commensurabilities it 
is possible for a portion of the vortices to be positioned 
at interstitial sites even for B/B^ < 1, so that the vor- 
tices form a type of ordering that would be absent for 
vortices in wire networks and Josephson-junction arrays. 
Further, the dynamic phases are also strongly affected by 
the fact that vortex flow can occur through the intersti- 
tial regions. 

For random pinning sites, theory [p4|-^, simulations 
|||,|^j28| and experiments (29|^ find that an initially 
pinned disordered vortex lattice can first flow plastically 
and then undergo a dynamic transition from a disor- 
dered structure to an ordered moving vortex lattice which 
can be a moving Bragg glass or moving smectic phase. 
The vortices in the strongly driven phase flow in well 
defined channels and a transverse depinning threshold 
can be present [p5| . Similarly, in simulations with peri- 
odic pinning where the vortex density exceeds the pin- 
ning density, B/B^ > 1, distinct kinds of plastic flow 



and large transverse barriers may exist 1 32 1 . Simulations 
of Josephson-junction arrays have also found plastic and 
elastic vortex flow phases and transverse depinning tran- 
sitions 11,11 . 

The dynamics of vortices in periodic pinning arrays 
with B/Bff, < 1 has not previously been studied. For 
this case, some of the open questions are (i) whether 
the vortex configurations in the pinned state are the 
same as the moving vortex configuration; (ii) how this 
depends on the filling fraction; and (Hi) whether there 
can be plastic fiow and how it depends on commensura- 
bility. Unlike systems with random pinning, where the 
disorder comes from the quenched substrate, the disor- 
der in vortex systems with periodic pinning comes about 
from the vortex-vortex interactions with the incommen- 
surate fields forming a disordered vortex lattice. At the 
commensurate fields, where the vortex lattice is ordered, 
plastic flow may be absent even though the vortex lat- 
tice is very strongly pinned. At strong driving forces it 
may be interesting to see if the incommensurate vortex 
lattice can reorder to a moving crystal phase in a similar 
manner as observed in systems with random pinning. It 
may also be interesting to compare commensurate and 
incommensurate fields in strongly driven phases to see 
whether the vortex flow occurs in ID channels and if the 
flow is between or along the pinning rows. It may fur- 
ther be valuable to correlate the microscopies of the static 
and dynamical phases to macroscopic observables, such 
as critical currents, current-voltage characteristics, and 
magnetization. 

In this work we examine the vortex configurations and 
dynamics for square and triangular pinning arrays in 2D 
superconductors with logarithmically interacting vortices 



using molecular dynamics simulations of pancake vor- 
tices. For the square pinning array we find a series of 
sub-matching effects that can be seen as peaks in the crit- 
ical depinning force. The vortex configurations at these 
fields agree well with those seen in recent imaging exper- 
iments [ p^p^Jl^ . Certain sub-matching fields, such as 
1/3 are not completely ordered, but are broken up into 
domains as also observed in experiments . Vortex or- 
dering at 1/3 filling depends on the pinning strength with 
stronger pinning producing a partially diagonally ordered 
arrangement. For weak disorder at 1/3 filling we observe 
a vortex ordering in which certain vortices sit in the in- 
terstitial regions between the pinning sites so that the 
overall vortex lattice has a near triangular ordering. We 
also find that the critical depinning force increases lin- 
early with pinning force strength at the commensurate 
fields while for low pinning strengths the incommensu- 
rate field shows a nonlinear (sub-linear) increase due to 
vortex-vortex interactions which crosses over to a linear 
increase. 

We find that at B/B^ = 1,1/2, and 1/4, where the 
vortex lattice is ordered, the depinning process and sub- 
sequent fiow is elastic with the moving vortex lattice hav- 
ing the same symmetry as the pinned lattice. For fields 
near strong commensurability, and for B/B^ > 1, we 
observe a two-stage depinning process where vacancies 
or interstitial defects in the commensurate vortex crystal 
depin first, as can be seen in velocity versus driving force 
curves that are equivalent to experimental V{I) curves. 

At the incommensurate fields, for B/B^ > 0.35, the 
vortex lattice depins plastically with only a portion of 
the vortices moving and the vortex motion occurring in 
winding paths in both the x and y directions. For these 
fields (vortex densities) at strong drives we observe re- 
ordering to an elastic distorted triangular moving lattice 
where the vortices move in ID channels along the pin- 
ning rows. For B / B^ < 0.35 the initial depinning is still 
plastic but is followed by a transition to a channel fiow 
phase where the vortices move in ID channels along the 
pinning rows. In contrast to what occurs at higher fields, 
the individual ID channels at lower fields contain differ- 
ent numbers of vortices, and the moving vortex lattice 
retains considerable positional disorder. These reordered 
phases are only stable in the strongly driven limit and 
we do not observe any hysteresis in the V{I) curves or 
vortex lattice order for either field range. 

For triangular pinning arrays we find a different set 
of matching filling fractions than those found for square 
pinning. The effective pinning of the lattice is enhanced 
and an ordered vortex lattice appears at B / B^ = 1, 6/7, 
2/3, 1/3, 1/4, 1/6, and 1/7. Further, there is no peak 
in the critical depinning force at 1/2. For strong pinning 
the vortex lattice at B/Bfp = 1/2 is disordered and all 
the vortices are located at pinning sites, while for weak 
disorder the vortex lattice at 1/2 has a square ordering 
in which half of the vortices are located at the inter- 
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stitial regions. For weak pinning the commensurability 
peaks in the critical depinning force are enhanced and an 
additional peak at B/B^ = 1/2 is observed. For both 
strong and weak pinning the moving vortex lattice at 
B/Bif, = 1/2 is square. For the triangular pinning array 
we do not find any matching effects for B/B^ > 1 which 
is due to the fact that the interactions with interstitial 
vortices cause a certain portion of the pinning sites to 
be left unoccupied. We observe V{I) characteristics and 
moving phases similar to the ones found in the square 
pinning array. 



II. SIMULATION 

We model logarithmically interacting pancake vortices 
in a 2D superconductor interacting with square and tri- 
angular pinning arrays. The overdamped equation of mo- 
tion for a vortex i is 



dt 



fVp 



fd = Vi 



The vortex-vortex interaction potential is Uy 



(1) 
ln(r). 

The total force force on vortex i from the other vortices 
is /™ = —Y^^^i^ iUvirij). We impose periodic bound- 
ary conditions and evaluate the periodic long-range log- 
arithmic interaction with an exact sum |3^. Pinning is 
modeled as attractive parabolic wells of radius rp, 



f7 = - Y.^U/r,){r. - r^') e((rp - |r, - r^^' |)/A) 



k=l 



(2) 



Here, we measure distance in units of the penetration 
depth A, is the location of pinning site fc, fp is the 
maximum pinning force, and Q is the Heaviside step func- 
tion The pinning is placed in a square or triangular array. 
The initial vortex positions are obtained by annealing 
from a high temperature configuration, with the vortices 
in a molten state, and slowly cooling to T = 0. This 
is intended to mimic field cooled experiments performed 
in We reduce the initial temperature to zero in 

twenty steps and remain at each step for 1.5 x 10^ MD 
steps, each being dt = 0.01 in normalized units. The vor- 
tex configurations are insensitive to the waiting time for 
waiting times greater than 5 x 10**. We note that the true 
ground states may only be attainable for prohibitively 
long waiting times; however, the agreement between our 
results and experiments is encouraging. 

After the vortex configurations are obtained, the criti- 
cal depinning force is determined by applying a slowly in- 
creasing, spatially uniform driving force, fd, which would 
correspond to a Lorentz force from an applied current. 
For each drive increment we measure the average vortex 
velocity in the direction of drive, Vx = Y^^" ^'^i- The 



force, /rf, versus velocity, T4, curve corresponds experi- 
mentally to a voltage-current, V{I), curve. The depin- 
ning force is defined to be the force value at which Vx 
reaches a value of 0.03 times the ohmic (linear) response. 



III. FRACTIONAL MATCHING FOR SQUARE 
PINNING ARRAY 

In Fig. l(a,b) we show the critical depinning force, /p, 
versus the vortex density for a system with a pinning 
density of tip — 1.36/A^, /p = 0.9/o, and Vp — 0.3A. In 
Fig. 1(a) we show /p for magnetic fields up to B/B^ = 
2.25, illustrating the matching effects at B/B^ — 1 and 
2 as well as clear matching effects at 1/2 and 1/4. Fig- 
ure 1(b) displays /^^ for B/B^ < 0.6, where additional 
sub-matching peaks are visible at 1/3, 1/6, 1/5, and 1/8, 
with the 1/3 matching field being the largest of these. 
We do not observe any particular matching at 2/3 or 4/5. 
Figure 2(a) shows the vortex positions (black circles) and 
pinning sites (open circles) for B / B^ = 1 , where the vor- 
tices form a square lattice with the pinning sites. In 
Fig. 2(b), we show the vortex configuration at an incom- 
mensurate field, B/B^ =^ 0.64, where a peak in the criti- 
cal depinning force is absent. Here, no apparent ordering 
in the vortex lattice can be observed. We also find that 
the vortex lattice is in general disordered at the other in- 
commensurate fields. In Fig. 2(c), for B/B^ — 1/2, the 
vortices form an ordered checkerboard state of a single 
domain. This state was also observed in imaging experi- 
ments on vortices in superconductors 1 10 l^Jl^ as well as 
of vortices in wire networks p6| , |37[ |. In the experiments 
by Field et al. ||l2| , domains of two different orientations 
of the checkerboard state are found at B/B^ = 1/2 where 
the domains of a single orientation are very large. It is 
unclear whether these domain structures in experiments 
result from some disorder at the pinning sites or are an 
intrinsic part of the vortex configurations as suggested by 
simulations of Josephson-junction arrays We have 

conducted larger simulations on systems with up to 392 
vortices at B/B^ — 1/2; however, we have found only 
single domains of vortices. In Fig. 2(d) the B/B^ = 1/3 
state is shown. Here the vortex configuration is not 
completely ordered; however, the vortices are generally 
positioned in diagonal stripes running along ±45° sepa- 
rated by two empty diagonal stripes. This state looks al- 
most identical to the results found in experiments ||l2|,|l3| . 
Completely ordered states at the 1/3 fiUing have been 
predicted for repulsive particles on a square lattice as well 
as vortices in square networks |^,^ . The disorder in our 
system may come from having too short of an annealing 
time. We have tested annealing times from 2.5 x 10^ to 
5 X 10^ normalized time units for B/B^ = 1/3 and have 
not found any appreciable differences in the vortex con- 
figurations; neither have we found noticeable system size 
effects. We find that for weaker pinning, the diagonal 
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ordering becomes less prominent as some vortices begin 
to sit in the interstitial regions suggesting that, for our 
simulations and possibly for the experiments, the vortex- 
vortex interaction, which tends to make the vortex lattice 
triangular, causes the disorder. In subsection B we show 
more explicitly the role of the pinning strength on the 
vortex lattice ordering at the 1/3 matching field. 

In Fig. 2(e) at B/B^ = 1/4 the vortices form an or- 
dered lattice with vortices occupying every other pinning 
site in every other row. In Fig. 2(f) we show the vor- 
tex lattice ordering at 1/5 which shows only regions of 
the square lattice ordering such as observed in the upper 
right area of Fig. 2(e). A similar ordering was observed 
in HI]. 

A. V(/) Characteristics 

We show, in Fig. 3(a), a series of V{I) curves for 
B/B^ = 0.98, 1, 1.05, and 1.2. The V{I) curve shows a 
single jump at the matching field. This is due to all the 
vortices depinning simultaneously. For B/B^ = 0.98 we 
observe a two-stage depinning process with initial depin- 
ning due to the onset of vacancy motion in the B jB^ = 1 
commensurate vortex lattice, followed by depinning of 
the vortices at the pinning sites which is associated with 
the second jump in the vortex velocities. AtS/B^ = 1.05 
the V{I) curves look similar to those obtained in previ- 
ous simulations for vortices in bulk superconductors . 
Here the pinned phase is considerably reduced compared 
to the 1/1 filling, and a two stage depinning process can 
be seen. The initial response is due to depinning of the 
interstitial vortices, and the second larger jump occurs 
when all the vortices begin to move. The additional 
jumps in the vortex velocities are due to transitions be- 
tween different plastic flow phases which are discussed in 
more detail in Ref. [|o). For B/B^ — 1.2 the two stage 
depinning features are still present and in general, for all 
B / B^ > 1, we observe the two depinning thresholds. 

In Fig. 3(b) we show velocity versus driving force for 
B/B^ = 0.525, 0.5, and 0.485. At 1/2 there is a sin- 
gle jump in Vx while at 0.525 and 0.485 the depinning 
transition occurs in two steps as seen in the initial jump 
in the V{I) followed by another larger jump. We find 
that two stage depinning occurs just above or below the 
matching fields where two distinct species of vortices ex- 
ist: the vortices forming the commensurate structure at 
the matching fields, and the vacancies or interstitial vor- 
tices in the commensurate vortex structure. These va- 
cancies and interstitials will have well defined depinning 
thresholds lower than those of the vortices forming the 
commensurate structure. For increasing or decreasing 
vortex density, relative to the 1/2 matching field, the 
overall vortex structure can change and the distinction 
between the two species of vortices is lost. In this case 
the two stage depinning transition disappears as seen in 



the V{I) curves in Fig. 3(c) for B/B^ = 0.3 and 0.41. 
Here the two-stage depinning process is clearly absent. 

B. Dependence of Vortex Configurations on Pinning 
Strength 

We have carried out a series of simulations at B/B^, = 
1, 1/2, 1/3, and 0.43 with fp ranging from O.l/o to 2.0/o. 
The vortex configurations at B/B^ = 1 and 1/2, as seen 
in Fig. l(a,b), are stable for all values of fp we have 
examined. In principle, even at these fillings, for small 
pinning strength, the elastic interactions of the vortex 
lattice will dominate, and consequently the lattice will 
be triangular. 

For B = 1/3 and for fp > O.6/0 the vortex configu- 
ration is the partially disordered phase with alternating 
diagonal rows filled as seen in Fig. 2(d). For weaker fp a 
portion of the vortices become located in the interstitial 
regions between the pinning sites. In Fig. 4(a) we show 
the vortices and the pinning sites and in Fig. 4(b) we 
show the vortices only for fp < 0.5/o where we observe 
an ordering when the vortex-vortex interaction is strong 
enough to induce triangular ordering; however, the vor- 
tex lattice can still take advantage of the periodicity of 
the pinning sites by allowing every other vortex to be 
pinned. In the ordered vortex lattice of Fig. 4(a), vor- 
tices fill every third pinning site in every other pinning 
row, while in the adjacent rows the vortices are not lo- 
cated in the pinning sites but instead are positioned in 
every third interstitial location. This produces a more 
triangular ordering in the vortex lattice which is easier 
to see in Fig. 4(b). As the pinning strength, fp, is re- 
duced for B/Bcf, = 0.43, we observe some vortices sitting 
at interstitial positions; however, we only observe regions 
of triangular ordering for the weakest pinning strength. 

In Fig. 5 we plot the critical depinning force, f^, versus 
pinning strength, fp, for B/ B^ = 1/2 and 0.43. Here, fp 
increases linearly with fp for the entire range of values 
considered, indicating that the vortices are in the single 
vortex pinning regime. Since the pinned vortex lattice 
is symmetric, this can be understood by realizing that 
the vortex-vortex interactions cancel and the depinning 
threshold will therefore be the pinning force strength, fp. 
For B/ Bfjj = 0.43 the critical depinning force increases in 
a non linear fashion for low fp and then becomes linear 
for high fp. At B/B^ = 0.43 the vortex lattice is disor- 
dered and the vortex-vortex interactions do not cancel, 
resulting in fp always being less than fp. The critical 
depinning force for B/B^ — 0.43 is always less than that 
of 1 /2 filling with the maximum difference occurring near 
fp = 0.3. 
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IV. MOVING VORTEX STATES FOR SQUARE 

PINNING ARRAY 

Figure 6 shows the vortex positions, pinning sites, 
and trajectories for the strongly driven vortex lattice for 
B/B^ = 1 [Fig. 6(a,b)] and 1/2 [Fig. 6(c,d)]. Here we 
show that the moving vortex structures are the same as 
the pinned vortex structures seen in Fig. 2. The depin- 
ning at 1/1, 1/2, and 1/4 filling fractions occurs elas- 
tically with vortices retaining the same neighbors they 
had at depinning. The vortices are also seen to flow in 
ID channels along the pinning sites. 

For incommensurate fields, B/B^ > 0.35, the initial 
depinning is plastic with vortices wandering in both the 
X and y directions. A transition to an ordered moving 
vortex lattice is found at higher drives. In Fig. 7 we 
show the vortex flow patterns for increasing driving force 
at B/B^ = 0.64. Just above depinning at fd = 0.805 
[Fig. 7(a)], the vortex flow is plastic with the vortices 
wandering in both the x and y directions. For stronger 
drives, fa/fo = 0.85 [Fig. 7(b)], the vortices begin to flow 
in ID channels along the pinning sites; however, there is 
still considerable hopping from one channel to the next. 
For fd/fo = 0.9 [Fig. 7(c)], the vortex lattice is almost 
completely reordered with only occasional vortex jumps 
between adjacent channels. For fd = 1.0 [Fig. 7(d)], the 
vortex lattice is completely reordered and the vortices 
flow in strict ID channels along the pinning sites. In 
Fig. 8 we show the vortex positions for fd/fo = 0.805 
[Fig. 8(a)], where the vortex lattice is clearly disordered, 
and fd/fo = 1-0 [Fig. 8(b)], where the vortex lattice has 
a defect-free distorted square ordering. The reordered 
vortex lattice can be frozen into the zero drive state by 
abruptly shutting off the driving; however, if the drive is 
slowly decreased the vortex lattice will re-enter the plas- 
tic flow phase and disorder so no hysteresis is observed 
in the V{I) curves. 

For B/B^ < 0.35 we find that the moving vortex lat- 
tice does not completely reorder. In Fig. 9(a) we show 
the vortex positions, pinning sites, and trajectories for 
B/B^ = 0.28 which indicate that the vortices are flow- 
ing elastically in ID channels along the pinning rows, but 
the number of vortices can differ from channel to chan- 
nel. In Fig. 9(b) the moving vortex lattice can be seen to 
retain a considerable amount of disorder. We have found 
that for weaker pinning the field at which the vortices 
can completely reorder is reduced below B/B^ = 0.35. 

V. FRACTIONAL VORTEX MATCHING AND 
DYNAMICS FOR TRIANGULAR PINNING 

ARRAYS 

In Fig. 10 we show the critical depinning force versus 
field for a system with triangular pinning. Matching ef- 
fects can be seen at B/B^ = 1, 6/7, 2/3, 1/3, 1/4, and 



1/7. No depinning peak is seen for the 1/2 matching nor 
ioY B/Bcj) > 1. In Fig. 11(a) we show the matching fields 
at B/B^ = 1, where all the vortices are pinned forming 
a triangular lattice. In Fig. 11(b), for B/B^ = 2/3, we 
find a honeycomb vortex lattice with vacancies forming a 
triangular lattice, which is the same as the vortex lattice 
configuration at the 1/3 matching field [Fig. 11(d)]. In 
Fig. 11(c), at B/B^ = 1/2, no particular vortex lattice 
ordering is observable. In Fig. 11(e), for B/B^ = 1/3, 
the vortices form a triangular sublattice. In Fig. 11(f), at 
B/B^ = 1/4, the vortices again form a triangular lattice 
while in Fig. 11(e), at B/B^ = 1/7, the overall vortex lat- 
tice is mostly triangular with some defects. The vortex 
lattice configuration at B/B^ = 1/7 (not shown) has a 
triangular configuration with a vacancy lattice the same 
as the vortex lattice in Fig. 11(e). 

A. Vortex Pinning and Ordering for Weaker Pinning 
Strength 

As was shown in Fig. 11(c), the vortex configuration 
at B/B^ = 1/2 is disordered. In Fig. 12(a) we show 
that, for a weaker pinning strength of /p//o = 0.3, an 
ordered state, in which half of the vortices are located 
in the interstitial regions between the pinning sites and 
the other half arc located at the pinning sites, is possible 
for B/B^ = 1/2. This gives the overall vortex lattice a 
square symmetry which is evident in Fig. 12(b). 

In Fig. 13 we show versus vortex density for fp/ fo = 
0.3 for a more limited field range ofO < B/B^ < 0.6. 
Here is considerably reduced from that of fp/fo = 0.9; 
however, the commcnsurability peaks are stronger with 
a larger difference between the commensurate peaks and 
the incommensurate regions. In addition to finding the 
same commcnsurability peaks observed at B/B^ = 1/3, 
1/4, and 1/7, we now clearly observe peaks at 1/6, 1/9, 
and 1/2. Most experiments in periodic pinning arrays 
have found that commcnsurability effects are strongest 
very near Tc where the pinning is the weakest. For de- 
creasing temperature, the overall pinning increases, but 
the clear matching anomalies gradually vanish. One pos- 
sibility (in the experiments) is that the bulk pinning be- 
comes relevant for lower T; however, experiments with- 
out dots or holes, for the same temperature regions, find 
that the critical current is much lower than in the system 
with the pinning arrays, indicating the strong pinning is 
coming from the dots or holes. Our simulations suggest 
that the increasing pinning force of the dots alone can 
wash out the matching anomalies even in the absence of 
bulk pinning. Figure 8 shows that the relative difference 
between the for a commensurate and incommensurate 
fields is maximum for low fp and is constant for higher 
fp. The relative height of the commensurate peak can be 
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expressed as 

^ (/p-(com) - /p-(inconi)) ^ 
^ /p(incom) /^(incom) ' 

For large fp, will be constant and /p(inconi) oc /p so 
5fp will decrease as ~ 1/ fp. 

The peak at 1/2 occurs due to the formation of the 
ordered vortex lattice as seen in Fig. 12. An experimen- 
tal test for the presence of this phase would be the ap- 
pearance of an anomaly at B/B^ = 1/2 in a triangular 
pinning array as T approaches T^. 

B. Vortex Ordering for 1 < B/B^y < 1.5 

In Fig. 14 we show the vortex configuration and pin- 
ning sites at B/B^ = 4/3, where it might be expected 
that the interstitial vortex lattice will resemble the vortex 
configuration at the B/B^ = 1/3 matching field. Fig. 14 
shows that a number of pinning sites are actually unoc- 
cupied with two vortices sitting adjacent to the pinning 
sites and caged in by six adjacent vortices. This indi- 
cates that placing an interstitial vortex between two vor- 
tices at the pinning sites is costly enough to depin one of 
the vortices. These types of defects disrupt any kind of 
overall interstitial vortex lattice ordering so that higher 
field fractional matching effects cannot be observed. For 
stronger pinning sites all the pinning sites will be occu- 
pied so fractional matching may be possible. 

The general V{I) characteristics are similar to those 
found in the square pinning arrays with a two step depin- 
ning process for B/B^ > 1 and within 0.05 flux-density 
of the sub-matching fields. We also find that the gen- 
eral features of the vortex dynamics and reordering are 
similar to those found with the square pinning array. 

VI. CONCLUSION 

We have investigated fractional matching effects in 
thin-film superconductors with square and triangular ar- 
rays of pinning sites. For the square arrays we find 
matching effects in the form of enhanced critical depin- 
ning forces at B/B^ = 1, 3/4, 1/2, 1/3, 1/4, and 1/5. 
The observed vortex configurations at these fields are in 
good agreement with recent imaging experiments. The 
vortex configurations also depend on the strength of the 
pinning sites. The vortex configuration at the 1/3 match- 
ing field can have a transition from a diagonally ordered 
lattice to a distorted triangular lattice, where every other 
vortex is interstitially pinned rather than being located 
in a pinning site. 

We also find that the critical depinning current scales 
linearly with fp. Since the vortex lattice at the commen- 
surate fields is symmetric, the vortex-vortex interactions 



cancel and the vortex lattice is always in the individual 
vortex pinning regime. For large fp, the critical depin- 
ning force for the incommensurate fields also scales lin- 
early with deviations from linearity at low fp as vortex- 
vortex interactions become relevant. 

Near B/B^ = 1,1/2, and 1/4 the V{I) characteristics 
show a two-stage depinning which is due to the presence 
of two well defined vortex species: the vortices of the 
commensurate lattice configurations, and the interstitials 
or vacancies in that lattice, which have different depin- 
ning thresholds. For B/B^ > 1, a two stage depinning 
in the V{I) curves is always observed. 

We find that the vortex lattice depins elastically at 
commensurate fields where the vortex lattice is symmet- 
ric; e.g., at B/B^ = 1, 1/2, and 1/4, while at the incom- 
mensurate fields the vortex lattice first undergoes plastic 
flow where the vortices flow in a random manner, and 
then at a higher drives reorders and flows elastically in 
well defined ID channels along the pinning sites. For 
B/B^ < 0.35, vortices will still flow elastically along the 
pinning sites; however, there is considerable disorder in 
the vortex lattice. For weaker pinning the vortex lattice 
will reorder for flelds B/B^ < 0.35. 

We have also investigated triangular pinning arrays 
and flnd matching effects at B/B^ = 1, 6/7, 2/3, 1/3, 
1/4, 1/5, and 1/7, which show ordered vortex lattices, 
while 1/2 and other incommensurate fillings show disor- 
dered vortex lattices. For B/B^ < 0.5 and for weaker 
pinning, we observe additional matching at 2/7 and 1/6 
as well as a general enhancement of the other commen- 
surate peaks. For the weaker pinning we also find a com- 
mensuration peak at B/B^ = 1/2, which occurs due to 
the formation of an ordering state where half the vortices 
are located in the interstitial regions so the overall vor- 
tex arrangement is square. The presence of these ordered 
states should be observable in magnetization curves. 

For triangular pinning arrays, with B/B^ > 1. we do 
not observe any matching features due to certain pinning 
sites remaining vacant where a pair of vortices can rotate 
around the vacant pinning site. We observe similar V{I) 
characteristics and dynamics for the triangular pinning 
array as seen in the square arrays. 

The static states for the pinning arrays should be ob- 
servable with Hall-probe arrays. Bitter decoration and 
Lorentz microscopy. Similar techniques may also be able 
to observe the dynamics; particularly at strong drives 
when the vortices should move in well defined ID chan- 
nels along the pinning. The motion of interstitials or 
vacancies near commensurate fields can be tested exper- 
imentally by looking for multiple depinning stages in the 
current-voltage curves. 
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FIG. 1. (a) Critical depinning force fp versus vortex den- 
sity B/B^ for < B/B^ < 2.1 with a square pinning ar- 
ray. Peaks can be seen in /p at B/B^ = 1/4, 1/2, 1/1, and 
2/1. (b) Critical depinning force for the system in (a) for 
< B/B^ < 0.6 showing narrowed view of the peaks in fp 
for B /B^ — 1/5, 1/4, 1/3, and 1/2. Smaller peaks can also be 
seen at B/B4, = 2/7, 1/6, and 1/8. 
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(e) (f) 
FIG. 2. Vortex positions (black circles) and pinning 

sites (open circles) from the simulation in Fig. 1 for (a) 

B/B^ = 1/1, (b), 0.642, (c) 1/2, (d) 1/3, (e) 1/4, and (f) 

1/5. 
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FIG. 3. Vortex velocities versus the applied driv- 
ing force fd for (a) B/B^ = 0.98,1.0,1.06, and 1.23. The 
largest critical depinning force occurs for the 1/1 matching, 
(b) B/B^ = 0.525, 0.5, and 0.485. (c) B/B4, = 0.30 and 0.41. 
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(a) (b) 
FIG. 4. (a) Vortex positions and pinning sites, (b) Vortex 
positions only for B/B^ = 1/3 with /p//o ~ 0.3. Unlike the 
case for for /p//o = 0.9 (Fig. 2(e)), where all the vortices are 
located at pinning sites, roughly half the vortices are located 
at the interstitial regions between the pinning sites so that 
the overall vortex structure is a distorted square as seen in 
(b). 




Pinning Strength fp/f„ 
FIG. 5. Critical depinning force versus fp for B/B^ = 1/2 
(upper curve) and B/B^ = 0.642 (lower curve). For 1/2 filling 
fp goes linearly with fp for the entire range of fp investigated. 
For 0.642 filling fp is linear in fp for large values of fp but 
is non-linear for low fp. A maximum difference between the 
two curves occurs near fp/fo = 0.4. 
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(c) (d) 
FIG. 6. (a) Vortex positions, pinning sites, and trajecto- 
ries for the moving vortex lattice at B/B^ = 1. (b) Vortex 
positions showing the moving square vortex lattice, (c) Vor- 
tex positions, pinning sites and trajectories for the moving 
system at B/B^ = 1/2. (d) Vortex positions showing the 
square moving vortex lattice rotated at 45° with respect to 
the pinning lattice. In both cases the moving lattice has the 
same symmetry as the pinned lattice. Further, the vortex 
motion can be seen to flow in well deflned ID channels along 
the pinning sites. 
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(a) 



(b) 




(c) 



(3) 

FIG. 7. Vortex positions, pinning sites, and trajectories 
for B/Bcp = 0.6428 for increasing applied driving force fd. 
(a) For fd/fo = 0.8, just above depinning, plcistic flow oc- 
curs with vortices flowing in both x and y directions, (b) At 
fd/fo = 0.9 the vortex motion is increasingly along the pin- 
ning sites forming ID channels, (c) At fd/fo ~ 0.95, almost 
all vortex flow is along the pinning sites with only occasional 
jumping of vortices between channels, (d) At fd/fo = 1, the 
vortex lattice has completely reordered. Vortices flow plasti- 
cally with the flow restricted to ID channels along the pinning 
sites. 



FIG. 8. Vortex positions from Fig. 7 for (a) fd/fo = 0.8 
and (b) fd/fo = 1-0. The strongly driven vortex lattice re- 
orders to a distorted square lattice. 
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FIG. 9. (a) Pinning sites, vortex positions and trajectories 
for B/B^ = 0.28 for fd/fo ~ 1. The vortices flow clastically 
in ID channels; however, the number of vortices is not equal 
in each channel, (b) Snapshot of the vortex positions, showing 
that the moving vortex lattice retains a considerable amount 
of disorder. 
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FIG. 10. Critical depinning force fp versus vortex density 
for a triangular pinning array. Commensuration peaks can be 
seen at B/B^ = 1/1, 6/7, 2/3, 1/3, 1/4, and 2/7. Note that 
there is no peak at B/B^ = 1/2. 
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(e) (f) 
FIG. 11. Vortex positions and pinning sites from Fig. 10 

for (a) BIB4, = 1/1, (b) 2/3, (c) 1/2, (d) 1/3, (e) 1/4, and 

(f) 1/7. 
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(a) (b) 
FIG. 12. (a) Vortex positions and pinning sites at 
B/B^ = 1/2 and (b) vortex positions only for fp/fo = 0.3. 
In (a), only half of the vortices are located at pinning sites 
with the other half located at the interstitial regions. In (b) 
the overall vortex lattice has a square symmetry. 
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FIG. 13. Critical depinning force for fp/fo = 0.3 for 
< B/B4, < 0.6. Commensurate peaks are observed 

at B/B^ = 1/4,1/3,1/6,2/7,1/7, and 1/9. A peak at 
B/B^ = 1/2 is now visible with the vortices forming the or- 
dered state observed in Fig. 10. 
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FIG. 14. Vortex positions and pinning sites at 

B/B^ = 1.25 showing that several of the pinning sites are 
unoccupied. 
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